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1. INTRODUCTION 

Consider the one-parameter family of integral equations of the second kind 

/•OO 

P(x,y) + ({x,y)+ dzP(x,z)ui(z,y) = 0, y > x, (1.1) 

J X 

where (3(x,y) is the unknown, ((x,y) is the nonhomogeneous term, and u(z,y) is an 
integral kernel which does not depend on the parameter x G R and satisfies 

/•OO 

sup / dz (\\u(z,y)\\ + \\u(y, z)\\) <+oo, (1.2) 

where || • || denotes any iV x N- matrix norm. Let us write (1.1) as 

/3 + C + /3O = 0, (1.3) 

where the integral operator Q acts from the right. From (1.2), as shown in the Appendix, 
it follows that O is bounded on the complex Banach spaces H^f xN of M x N matrix- valued 
measurable functions F : (x, +oo) — > C MxN such that the matrix norm ||F(-)|| belongs to 
L p (x, +oo) for 1 < p < +oo. 

We assume that, for each x G R, (7 + 0) is an invertible operator on H^ xN and on 
H.2 XN , where J denotes the identity operator. Using (I + R) to denote the corresponding 
resolvent operator, where 

R:= (I + O) -1 — I, 7 + i?= (I + fi)- 1 , (1.4) 

the solution to (1.3) can formally be written as 

P=-((I + R), 

or equivalent ly as 

/•OO 

P(x,y) = ~({x,y) - / dzC(s,z)r(i;z,|/), (1.5) 

where r(x; 2, y) denotes the integral kernel of the operator R. 
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Let us consider (1.3) in the special case 

a + w + aO = 0, (1.6) 

where the nonhomogeneous term and the integral kernel coincide, as seen by writing (1.6) 
explicitly as 

/•OO 

a(x, y) + u{x, y) + / dz a(x, z) u>(z, y) = 0, y > x. (1.7) 

J X 

The solution to (1.6) can formally be written as 

a = + R). (1.8) 

The unique solvability of (1.7) in 7i^ xN and the condition in (1.2) imply that 

/•oo 

sup / dz{\\a(z,y)\\ + \\ct(y,z)\\) < +oo. (1.9) 
y>x Jx 

A fundamental question related to (1.3) is the following: Can we write r(x;y,z) ap- 
pearing in (1.5) explicitly in terms of a(x,y) appearing in (1.7)? In case the answer is 
affirmative, we can express the solution (3 to (1.3) explicitly in terms of a and (. In fact, 
such a reduction question dates back to the Armenian astrophysicist Ambarzumian whose 
invariance principles are used in transfer of light in planetary atmospheres [7,11,12,32,33]. 
Ambarzumian [6] considered (1.6) with uj{y,z) = (c/2)Ei(|y — z\), where Ei is the ex- 
ponential integral function. Similar reduction formulas were obtained [20,21] for integral 
equations with convolution kernels, i.e. when cu(y, z) is a function of (y — z). 

One of our goals in this paper is to study the aforementioned fundamental question 
when the integral operator fi in (1.3) is N x A^-matrix valued and J-selfadjoint in the sense 
that 

n = J&J, u(y,z) = Juo(z,y) ] J, (1.10) 

where the dagger denotes the matrix adjoint (complex conjugation and matrix transpose) 
and J is an A x A selfadjoint involution, i.e. 



We present one of our key results in Theorem 2.2, where the resolvent kernel r(x;y,z) 
appearing in (1.5) is explicitly expressed in terms of the solution a(x,y) to (1.6). 



J : = 



Let us note that, without loss of generality, J may be assumed to have the form 



-In-j 

where Ij is the j x j identity matrix for some 1 < j < N. In that case we have 
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for block matrices Mi, M 2 , M 3 , M 4 of appropriate sizes. 

Having established our first key result in Theorem 2.2, we turn our attention to the 
integral equation 

a + cD + aO = 0, (1.11) 

which is explicitly written as 

pOC 

ct(x,y) + &{x,y) + / dz a(x, z) u>(z, y) = 0, y > x, (1-12) 

J X 

obtained from (1.6) by perturbing the operator O to O by a finite-rank operator, i.e. 

= + FG, u(x, y) = u(x, y) + f(x) g(y), (1.13) 

where / and # are Nxj and j xN matrices with entries depending on a single independent 
variable and belonging to 7if rxj fl 7i^ Xj and Tt{ xN fl W^; N ', respectively. We note that 
we cannot in general expect F and G to commute, and hence in general /<7 ^ In our 
second key result, we show that (1.11) can be transformed into another integral equation in 
which the kernel is degenerate (i.e. separable in the independent variables) so that a(x, y) 
can be explicitly obtained in terms of a(x, y), f(x), and g(y), as indicated in Theorem 3.4. 

Our key result given in Theorem 3.4 has important implications for various linear 
differential equations or systems of differential equations arising in important physical 
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applications. One important consequence of Theorem 3.4 is that it provides a systematic 
method to derive the Darboux transformations for a wide variety of spectral problems 
for differential equations. Recall that the idea behind a Darboux transformation (see e.g. 
[9,27,31] and the references therein) is to determine how the (generalized) eigenvectors 
change when a finite number of discrete eigenvalues are added to or subtracted from the 
spectrum of a differential operator without changing the continuous spectrum. In the 
language of physics, the Darboux transformation provides the perturbed potential and 
wave function in terms of the unperturbed quantities when a finite number of bound 
states are added or subtracted. 

In this paper we are only concerned with Darboux transformations and not with 
Backlund transformations. When a discrete eigenvalue is added to the spectrum of a dif- 
ferential operator, a Backlund transformation [18] usually consists of a first-order differen- 
tial equation (or a system of first-order differential equations) involving the perturbed and 
unperturbed potentials. On the other hand, in a Darboux transformation the perturbed 
potential is explicitly expressed in terms of unperturbed quantities. Backlund transfor- 
mations have been derived [8,13,19,25,34] for various systems of differential operators and 
they are useful in obtaining exact solutions to related nonlinear evolution equations such 
as the Korteweg-de Vries equation, the nonlinear Schrodinger equation, the sine-Gordon 
equation, and the modified Korteweg-de Vries equation. 

The Darboux transformation is well understood for Sturm-Liouville problems on a 
finite interval [14] and the one-dimensional Schrodinger equation [15], but there are also 
many others where such transformations are not yet known or only some very special cases 
are known. Our Theorem 3.4 can be applied to the Zakharov-Shabat differential operator, 
matrix Zakharov-Shabat systems, and other differential operators to derive in a systematic 
way the corresponding Darboux transformations both at the potential and wave function 
levels. Our theorem is general enough so that it applies when one eigenvalue is added or 
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subtracted from the spectrum, several eigenvalues are added or subtracted simultaneously, 
and eigenvalues with nontrivial Jordan structures are added or subtracted either one at a 
time or simultaneously. As an example, we apply Theorem 3.4 on the Zakharov-Shabat 
system, and in Theorem 6.4 we present the Darboux transformation expressing both the 
change in the potential and the change in the wave function explicitly in terms of the wave 
function of the unperturbed problem when one bound state is added. We compare our 
transformation for the potential given in (6.18) and for the wave function given in (6.19) 
with the results in the literature. 

Let us mention that our results remain valid if the range of the integral is over (— oo, x) 
so that (1.7) is replaced with 

a(x, y) + u(x, y) + / dza(x, z) u>(z, y) = 0, y < x, (1-14) 

J — oo 

and also valid if the integral is over (0, x) so that (1.7) is replaced with 

a(x,y) + u(x,y) + / dz a(x, z) u)(z, y) = 0, < y < x, (1-15) 
Jo 

and with the obvious appropriate replacements in (1.1), (1.2), (1.5), (1.7), (1.9), and (1.12). 
By using the operator notation of (1.3) and (1.6) it is straightforward to modify the proofs 
and to treat (1.7), (1.14), and (1.15) all at once. 

Our paper is organized as follows. In Section 2 we establish our first key result 
by expressing the resolvent kernel r(x;y,z) appearing in (1.5) explicitly in terms of the 
solution a(x, y) to (1.6). In Section 3 we obtain our second key result by expressing a(x, y) 
explicitly in terms of a(x,y), f{x), g(y) when (1.6) is perturbed to (1.11) as in (1.13). In 
Section 4 we show how the key result in Theorem 3.4 provides a unified approach to 
derive Darboux transformations. In Section 5 we show that the key results in the previous 
sections are applicable to various systems such as the Zakharov-Shabat system, its matrix 
generalizations, and the Schrodinger equations on the full and half lines; we show how an 
integral equation of the form (1.7), (1.14), or (1.15) arises for each system and is related 
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to an associated Marchenko integral equation or a GePfand-Levitan integral equation. 
In Section 6, we illustrate the significance of our Theorem 3.4 and derive the Darboux 
transformation for the Zakharov-Shabat system and make a comparison with some related 
results in the literature. 

2. REDUCTION OF THE RESOLVENT KERNEL 

Recall that we assume that (1.3) is uniquely solvable in Hi XN and in 7i^ xAr and that 
the operator O and its integral kernel u(y,z) satisfy (1.2) and (1.10). In this section, we 
analyze the resolvent kernel r(x;y,z) appearing in (1.5) and present our first key result; 
namely, we show that r(x; y, z) can be expressed explicitly in terms of the solution a(x, y) 
to (1.7). 

Proposition 2.1 Assume that (1.3) is uniquely solvable in 7i^ xAr and that O satisfies 
(1.10). Then, the operator R given in (1-4) and the corresponding kernel r(x;y, z) appear- 
ing in (1.5) satisfy 

R = JR i J, r(x;y,z) = Jr{x;z,y)U, (2.1) 
where J is the involution matrix appearing in (1.10). 
PROOF: From (1.4) we see that 

(J + Q)(I + R) = / = (/ + R)(I + O), 

and hence we obtain 

R + n + QR = 0, (2.2) 

R + Q + RQ = 0. (2.3) 

By taking the adjoint of the operator equation in (2.2) and applying J on both sides of 
the resulting equation, we get 

JR^J + JO t J + (JR)J){JQ}J) = 0, 
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or equivalently, after using (1.10), 

JR^J + Q + {JB)j)Vt = 0. (2.4) 

Since (1.3) is assumed to be uniquely solvable in 7Y^ xAr , by comparing (2.3) and (2.4) we 
see that R = JR^ J. In taking the adjoint, we note that the independent variables y and z 
are switched in the argument of the kernel and hence (2.1) is established. | 

Our first key result is given in the next theorem. 

Theorem 2.2 Assume that (1.3) is uniquely solvable in 7i^ x7V and that O satisfies (1.10). 
Then, the corresponding kernel r(x;y, z) appearing in (1.5) can be expressed explicitly in 
terms of the solution a(x,y) to (1.7) as 

a(y,z)+ / ds J a(s, yy J a(s, z), x < y < z, 
r(x;y,z)=l ' ^ (2.5) 

Ja(z,y)^J+ / ds J a(s,yy J a(s, z), x < z < y, 

V Jx 

where J is the involution matrix appearing in (1.10). 

PROOF: Since (1.3) is uniquely solvable in 7i^ x7V , so is (1.6) and hence the solution R to 
(2.3) is unique. Thus, it suffices to prove that the quantity defined in (2.5) satisfies (2.3), 
i.e. the quantity in (2.5) satisfies the integral equation 

/•OO 

r(x;y, z) + oo{y, z) + / ds r(x; y, s) u(s, z) = 0, a; < min{y, z}. (2.6) 

J X 

The proof for the case x < z < y is similar to the case x < y < z, and hence we will only 
give the proof in the latter case. In that case, let us use = + in the integral 
appearing in (2.6). We use a direct substitution from (2.5) into (2.6), where we note that 
the first line of (2.5) is used in the integral and the second line of (2.7) is used in the 
integral in (2.6). After the substitution, the left hand side in (2.6) becomes i>i +i>2 + ^3, 
where we have defined 

/•OO 

v 1 :=a(y,z) + u(y,z)+ / ds a(y, s) u(s, z), 

Jy 



ry ry 
V2 := / dt J a(t, yy J uj{t, z) + / dt J a(t,yy J a(t, z), 

J X J X 

ry rs roc py 

V3 := ds dt J a(t,y)^ J a(t, s) oj(s, z) + / ds dt J a(t,y)^ J a(t, s) oj(s, z). 

J x J x J y J x 

Note that vi = because of (1.7). The orders of the two iterated integrals in v 3 can be 
changed to f£ dt ds and dt ds, respectively. Using + = J t °°, we then get 

ry r r°o 

V2 + vs = / dtJu(t,y)^J a(t, z) + oj(t, z) + / ds a(t, s) oj(s, z) . (2.7) 

Jx L Jt 

We see that the quantity in the brackets in (2.7) vanishes because of (1.7). Thus, (2.6) is 
satisfied for x < y < z. A similar direct substitution for the case x < z < y completes the 
proof. I 

3. FINITE-RANK PERTURBATIONS 

Our main goal in this section is to show that the solution a(x, y) to (1.12) can be 
expressed explicitly in terms of a(x,y), f(x), and g(y) appearing in (1.7) and (1.13), 
respectively. As we will see in later sections, the key formulas given in (3.19)-(3.21) below 
form the basis of Darboux transformations related to a wide variety of spectral problems 
associated with ordinary differential operators. 

Recall that we assume that (1.10) holds and that (1.7) is uniquely solvable on H^ xN 
and onH^ xN . Let us now define the intermediate quantities n(x) and q(y) as 

POO POO 

n(x) := f(x) + / dz a(x, z) f(z), q(x) := g(x) + / dz g(z) J a(x, z)^ J, (3.1) 

J X J X 

where J is the involution matrix appearing in (1.10). From (1.9) it follows that n G 
Hi xj D H^ xj and q G H{ xN fl H£ N . Note that both integration limits J x °° in (3.1) can 
be replaced with because a{x, y) = for x > y. 

Theorem 3.1 We can transform (1.11) into an integral equation that has a degenerate 
kernel and hence obtain a explicitly by linear algebraic methods. 
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PROOF: Using (1.13) let us write (1.11) as 

a(I + tt + FG) = -u - fg. (3.2) 

Recall that all the operators act from the right. By applying on (3.2) from the right with 
the resolvent operator (I + R) appearing in (1.4), we get 

a[I + FG(I + R)]=a-fg(I + R), (3.3) 

where we have used (1.8) to have a on the right hand side of (3.3). Let us define the 
operator G as 

/>oo 

G:=G(I + R), ~g{x,y):=g{y)+ dz g(z) r(x; z,y), (3.4) 

J X 

where r(x;y, z) is the kernel given in (2.5). We emphasize the dependence of g both on x 
and y. Note that the integral equation in (3.3) has a degenerate kernel, which can be seen 
by writing it in the form 

a(I + FG) = a — fg, (3.5) 

because the kernel of FG is f(y) g(x, z), where there is a separation of the y and z variables, 
and x appears merely as a parameter. 

Let us now solve (3.3) by using linear algebra. We look for a solution in the form 

a(x,y) = a(x,y) +p(x)g(x,y), (3.6) 
where p is to be determined. Using (3.6) in (3.5), after some simplification we get 

(aF+p+pgF + f)G = 0, 

which yields 

p(I+~ g F) = -(f + aF), (3.7) 
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or written in the integral form as 



p(x) = —n(x) 



poo 

1+ dsg(x,s)f(s) 

J X 



where we have used the definition of n(x) given in (3.1). Using (3.7) in (3.6) we obtain 



a = a-(f + aF)(I + gF)- l g 



or written in the integral form as 



a(x,y) = a{x,y) - n(x) 



1+ dsg(x,s) f(s) 

J X 



(3.8) 



which completes the proof of our theorem. | 



Note that (3.8) expresses a(x, y) in terms of a(x, y), /(x), and g(x, y) because as seen 
from (3.1) the quantity n(x) is available in terms of a(x,y) and f(x). 

Next, we show that g(x,y) can explicitly be obtained in terms of a(x,y) and g(y), 
which will then imply that a(x,y) is expressed in terms of a(x,y), f(x), and g(y). 

Proposition 3.2 The quantity g(x,y) defined in (3.4) can be expressed explicitly in terms 
of the solution a(x,y) to (1.6) and the quantities f(x) and g(y) appearing in (1.13) as 



g(x,y)=q(y)+ dsq(s)a(s,y), (3.9) 

J X 

where q(y) is the quantity defined in (3.1), and it is noted that g(x,x) = q(x). 



PROOF: We will substitute (2.5) into (3.4). For this purpose, let us write the integral J a 



oc 
x 



in (3.4) as + . Using the first line of (2.5) in and the second in f^° 7 we obtain 



9( x ^y)=9(y)+ dsg(s)a(s,y)+ ds g(s) J a(y, s) f J 

J x J y 

/ rV rs roc, py \ 

+ ( ds dt+ ds dt\ g(s) J 'a(t,s) f J a(t,y). 

\J x J x J y J x J 

The sum of the two iterated integrals above can be written first as a double integral and 



then as an iterated integral by changing the order of integration to get 



ry r s poo ry ry poo 

ds dt+ ds dt= dt ds. 

Jx Jx Jy Jx Jx Jt 



(3.10) 
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Using (3.10) and combining terms as in (3.1), we then obtain (3.9). | 

We note that the integral in (3.9) can also be written as because a(s,y) = 
for s > y. 

Let us define the matrix T(x) as the quantity whose inverse appearing in (3.8), namely 



as 



pOO 

T(x) :=/+ / ds~g(x,s)f(s). (3.11) 

J X 

Proposition 3.3 The quantity T(x) defined in (3.11) can be expressed explicitly in terms 
of the solution a(x,y) to (1.6) and the quantities f(x) and g(y) appearing in (1.13) as 

POO 

T(x)=I+ dsq(s)n(s), (3.12) 

J X 

where n(x) and q(x) are the quantities defined in (3.1). 
PROOF: Using (3.9) in (3.11) we get 

poo poo ps 

T(x) = I+ dsq(s)f(s)+ ds / dt q(t) a(t, s) f(s). (3.13) 

Jx Jx Jx 

Changing the order of integration in the last integral in (3.13) and using (3.11), we get 
(3.12). Since n(s) and q(s) are expressed in terms of a(x, y), f(x), and g(y), we see from 
(3.12) that T(x) is explicitly expressed in terms of a(x,y), f(x), and g(y), as well. | 

The Fourier transform of the N x N matrix-valued quantity a(x, y) in (1.7), usually 
called a wave function, can be written as 

/•OO 

*(A,x) :=e- %XJx + / dya(x,y)e- lXJ \ (3.14) 

J X 

where J is the involution matrix appearing in (1.10). Using the inverse Fourier transform 
on (3.14) we get 

1 f°° 

a (x, y) = 7r d\ [tf (A, x) - e ~ lXJx ] e tXJy . (3.15) 
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Similarly, for a(x, y), we have the associated N x N matrix-valued wave function ^(A, x), 
where 

/•OO 

V(\,x) :=e- lXJx + / dya(x,y)e- lXJy , (3.16) 

J X 



i Z" 00 r~ 



-iXJx 



Let us introduce 



poo 

7(A,x):= / dyS(x,y)e- <AJ ». (3.17) 
Using (3.9) and (3.14) in (3.17) and the fact that a(x, y) = for y < x, we get 

p OO 

7(A,x) = / dsg(s)*(A,s). (3.18) 

The following theorem and in particular (3.20) describes the effect of a finite-rank 
perturbation on the wave function. 

Theorem 3.4 Let a and a be the solutions to the integral equations (1.6) and (1.11), 

respectively, and let n(x), T(x), g(x,y), and^y(X,x) be the quantities given in (3.1), (3.12), 
(3.9), and (3.18), respectively. Then, a(x,y) — a(x,y) and ^(A, x) — ^(A, x) can explicitly 
be written in terms ofa(x,y), f(x), g(y) as 

&(x, y) - a(x, y) = -n(x) T(x)~ 1 g(x, y), (3.19) 

*(A, x) - *(A, x) = -7i(x)r(x)" 1 7(A, x). (3.20) 

Furthermore, we have 

a(x, x) — a(x, x) = —n(x) T(x)~ 1 q(x). (3-21) 



PROOF: Note that using (3.1), (3.11), and (3.12) in (3.8) we obtain (3.19). Using g(x, x) = 
q(x) from Proposition 3.2 in (3.19) we get (3.21). Finally, we obtain (3.20) with the help 
of (3.14), (3.16), (3.17), and (3.19). | 
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We conclude this section with a result on the trace of the left hand side of (3.21). 

Proposition 3.5 Let a and a be the solutions of the integral equations (1.6) and (1.11), 
respectively, and let Y be the matrix given in (3.12). The trace of the difference a(x,x) — 
a(x,x) can be expressed as the logarithmic derivative of the determinant ofT(x) as 



tr [a(x, x) — a(x, x)] = tr 



dT{x) 
dx 



T(x)- 1 



1 ddetT(x) 
detr(x) dx 



PROOF: From (3.12) we get 



dT(x) 
dx 



—q(x) n(x). 



Using the well-known matrix properties 



tr [MiM 2 ] = tr [M 2 Mi] , tr 



dM{x) 
dx 



M(x)- 1 



1 ddetM(x) 
det M(x) dx 



from (3.21) and (3.23) we get (3.22). | 



(3.22) 



(3.23) 



4. DARBOUX TRANSFORMATIONS 

The integral equations (1.7), (1.14), and (1.15) arise in the study of various scattering 
and spectral problems, some of which are described in Section 5. In this section we will 
elaborate on (3.19)-(3.21) and show how they provide a unified approach to Darboux 
transformations for a variety of scattering and spectral problems. Recall that a Darboux 
transformation describes how the wave function and the potential change when a finite 
number of discrete eigenvalues are added (or subtracted) from the spectrum of a differential 
operator without changing its continuous spectrum. 

Suppose we add a discrete eigenvalue Xj with multiplicity nj to the existing spectrum. 
Then, associated with the eigenvalue Xj, there are nj parameters cjq, . . . , Cj( n i), usually 
known as norming constants. The formulas (3.19)-(3.21) tell us how the wave function 
changes from \I/(A, x) to ^(A, x), how the potential changes from u(x) to u(x), and how the 
quantity a(x,y) related to the Fourier transform of the wave function changes to a(x,y). 
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Consequently, for each discrete eigenvalue Xj added to the spectrum, there will be an rij- 
parameter family of potentials u(x), where the norming constants act as the parameters. 
In case several discrete eigenvalues Ai,...,Ajv are added all at once, it is convenient to 
use [4] a square matrix A whose eigenvalues are related to Xj for j = 1, . . . , N in a simple 
manner; it is also convenient to use [4] a matrix C whose entries are related to the norming 
constants c JS for j = 1, . . . , N and s = 0, 1, . . . , rij — 1. 

The quantities f(x) and g(x) appearing in (1.13) can usually be represented in the 

form 



/(*) = 



Bh~ Atx 
Ce~ Ax 



, g(y) = 



e~ Ay B 
-e 







(4.1) 



where A is a constant square matrix with all eigenvalues having positive real parts (the 
bound-state A-values are usually obtained [4] by multiplying the eigenvalues of A by the 
imaginary unit i), and B and C are constant matrices of appropriate sizes so that the 
matrix product f(x) g(y) is well defined and given by 

-Bte-^aH-yJct' 

Ce -A(x+y) B q 



f(x)g(y) 



For f(x) and g(y) given in (4.1), let us evaluate n(x), q(x), and g(x,y) given in (3.1) 
and (3.9), respectively, explicitly in terms of the wave function ^(A, x) evaluated at the 
eigenvalues of A. First, by taking the matrix adjoint, from (3.15) we get 



1 f°° 

Ja(x, y)t J = — / dX e~ lXJy [J*(A, x) t J - e lXJx ] . 

^ n J — oo 

Using (3.15) in (3.1) and the fact that 

1 r°° 

— J_ dse^ as = 5(a), 



(4.2) 



where 5 is the Dirac delta distribution, we evaluate n(x) as 



1 poo POO 

n{x)= 2^J dX *^ x ^J dye* XJ yf(y). 



(4.3) 



(4.4) 
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Using (4.1) in (4.4) we obtain 

1 f°° 

n(x) = — / dXV(X,x) e iXJx N{X,x), 

where we have defined 

-BUM + iA^)- 1 e- A ' x 

_C(XI -iA)- x e- Ax 

Similarly, using (4.2) and (4.3) in (3.1) we obtain 

-I poo poo 

q(x) = — dX dyg(y)e-* XJ yjV(X,x)U, 

J-oo Jx 

and using (4.1) in (4.6) we conclude that 

1 f°° 

q( x ) = / dX Q(X, x)e~ lXJx J^(X, x) ] J, 

2?" J-oo 



(4.5) 



(4.6) 



(4.7) 



where we have defined 



Q(X,x) := 



e~ Ax (XI -iA^B 







e~ Atx (XI + iA^)~ 1 C\ 

Proceeding in a similar manner, with the help of (3.9), (3.14), (4.3), and (4.7) we first 



get 



^ poo poo 

9( X ^) = ^J d vJ dzq{z)^(ii,z)e^ J \ 



and then obtain 



/OO pOO poo 

dX dfi dzE(X,fi 7 z)- 
-oo J — oo Jx 



A^Jy 



(4.8) 



where we have defined 



E(X, (i, x) := Q(A, x) e~ lXJx Jtf (A, x) ] J Vfa, x). 



(4.9) 



Then, with the help of (4.5) and (4.7) the matrix T(x) given in (3.12) can be explicitly 
written in terms of the wave function \P as 

/oo poo poo 

dX d/i dyE(\,n,y) e^ Jy U{^y). (4.10) 
-oo J — oo Jx 
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Finally, using (4.7) and (4.9) in (3.18), we obtain 

/•oo />oo 

7(A,x) = 2tt/ ds d\E(X, A, s). (4.11) 

J x J — oo 

Note that the integrals in (4.5), (4.7), (4.10), and (4.11) can be performed as residue 
integrals in the complex A-plane with the poles at the eigenvalues of iA and —iA^. Evaluat- 
ing those integrals and using the result in (4.8) and (4.11) we can obtain g(x, y) and 7(A, x) 
explicitly in terms of \1/(A, x) evaluated at the eigenvalues of iA and —iA^. If some bound 
states have multiplicities, i.e. if some of the eigenvalues of A have nontrivial Jordan struc- 
tures, then the explicit expressions for g(x, y) and 7(A, x) also contain some A-derivatives 
of \1/(A, x) evaluated at the eigenvalues of iA and —iA^. 

In Section 6 we will use the procedure described here to obtain the Darboux transfor- 
mation for the Zakharov-Shabat system given in (5.1). 

5. APPLICATIONS TO SPECIFIC SYSTEMS 

In this section we present some specific systems on which the theory presented in the 
previous sections is applicable. In the first four examples we analyze the Zakharov-Shabat 
system and its matrix generalizations. In the remaining three examples we analyze the 
Schrodinger equation on the full and half lines. For each system, we identify the quantities 
uj(x,y) and a(x,y) appearing in (1.7) or in one of its variants (1.14) and (1.15), and 
we identify f(x) and g(y) appearing in (1.13). We identify the involution J for which 
(1.10) is satisfied by ou(y,z) for each system. Then, for each system we indicate how the 
potential u(x) appearing in the system is related to a(x, x) so that one sees clearly how the 
perturbation u(x) — u(x) in the potential can be recovered from (3.21). For each system 
we indicate how the wave function \I/(A, x) is related to certain specific solutions to the 
corresponding system. We also relate the integral equations (1.7), (1.14), and (1.15) to 
the Marchenko equations or the Gel'fand-Levitan equation corresponding to each system. 
In some cases we observe that (1.7), (1.14), or (1.15) is exactly the same as a Marchenko 
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equation or a Gel'fand-Levitan equation, and in some other cases one needs to rearrange 
the Marchenko equations or the Gel'fand-Levitan equation in order to get the integral 
equation (1.7), (1.14), or (1.15). 

We first discuss four examples involving the Zakharov-Shabat system and its matrix 
generalization, two with the range of the integral over (x, +oo) and two with the range of 
the integral over (— oo, x). We shall discuss the details of Example 5.1 in the next section, 
where we derive the Darboux transformation for the Zakharov-Shabat system and compare 
our results with the existing results in the literature [10,19,23,24,26,29]. 

Example 5.1 Consider the Zakharov-Shabat system 



d(p(X, x) 
dx 



—iX u(x) 
-u(x)* iX 



(f(X,x), 



(5.1) 



where the asterisk denotes complex conjugation, u is a (scalar) complex-valued integrable 
potential, and <p is a column vector with two components. The corresponding (left) 
Marchenko equations are given by [30] 



K{x,y) + 



K(x,y) 





Q\(x + y)_ 

ni(x + y y 





+ / dzK(x,z)Qi(z + y) = 



/•OO 

/ dzK(x,z)Qi(z + y) ] 

J X 



y > x, 



y > x, 



(5.2) 



where K(x, y) and K(x, y) are 2x1 matrix valued and Q\(x+y) is a scalar function. Let us 
stress that an over line does not indicate complex conjugation. The potential is recovered 
as [30] 



u(x) = -2[l 0}K(x,x) = 2K(x,x) ] 



(5.3) 



By letting 



a(x,y) = [K(x,y) K{x,y)}, u(x,y) = 



-fii(x + y)t 

Cli(x + y) 



(5.4) 



18 



we can write (5.2) as (1.7), which is now a 2 x 2 system of integral equations. We note 
that u{x,y) given in (5.4) satisfies (1.10) with J = diagfl, — 1]. As seen from (5.3), u is 
then recovered from the solution to (1.7) as 



u(x) = -2[l 0]a(x,x) 



(5.5) 



In this case, the degenerate perturbation on Q\(x + y) is given [4] by Cie- Ai(x+y ^B h where 
A\ is a constant pxp matrix with all eigenvalues having positive real parts, B\ is a constant 
p x 1 matrix, and C\ is a constant 1 x p matrix. The functions / and g appearing in (1.13) 
are then 2 x 2p and 2p x 2 matrices, respectively, given by 



B\e~ A \ l 
Cie~ AlX 



g(y) 



e~ Aiy B l 






-e- A \yc\ 



(5.6) 



For the Zakharov-Shabat system the wave function appearing in (3.14) is the 2x2 matrix 
given by 



*(A,x) = 



ip 2 (\*,x)* ipi(X,x) 
-V»i(A*,x)* V2(A,x) 



(5.7) 



where 
and 



^i(A, x) 
^ 2 (A,x 



is the Jost solution to (5.1) with the asymptotics 

e -i. 







as x — >■ +oo 



is the solution behaving as 







+ o(l) as x — > +oo. 



Example 5.2 Consider the matrix generalization of Example 5.1 with the Zakharov- 
Shabat system 

d<p(\,x) \-iUm "!•'•) I x 

¥>(A,x), (5.8) 



-iXIm U(x) 
—u(xY i\I n 

where u is an m x n matrix with integrable entries, and Ij is the j x j identity matrix. 
The corresponding (left) Marchenko equations are given by 



K(x,y) + 



K{x,y) - 



Omn 

Q\(x + y) ] 

rt\(x+ y y 

77.777. 



poo 



dzK{x,z) Qi(z + y) ] 



r 



nn 
On 
0„ 



2/ > x, 



y> x, 



(5.9) 
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where Ojk is the zero matrix of size j x k, and K, K, and Vt\ have sizes (to + n) x n, 
(m + n) x to, and n x to, respectively. The to x n potential matrix u is recovered from the 
solution to (5.9) as 

Orn/i 



(5.10) 



u(x) = -2[I m mn ] K(x, x) = 2K(x,x) ] 

By letting 

0mm -O i (x + y) t 

a(x,y) = [K(x,y) K(x,y)}, iv(x,y) = 

_lh(x + y) nn 

we can write (5.9) as (1.7), which is now an (m + n) x (m + n) system of integral equations. 
Note that ou(x, y) given in (5.11) satisfies (1.10) with J = I m @ (—In)- As seen from (5.10), 
the potential u is recovered from the solution to (1.7) as 



(5.11) 



u(x) 



-2[I m 



m7 



In this case, the degenerate perturbation on the n x to matrix quantity Q\(x + y) is given 
[4,16,17] by C\e~ Al ^ x+yS) B\, where A\ is a constant pxp matrix with all eigenvalues having 
positive real parts, B\ is a constant p x to matrix, and Ci is a constant n x p matrix. The 
functions / and g appearing in (1.13) are then (m + n) x 2p and 2p x (to + n) matrices, 
respectively, given by 



fix) 







mp 



Cie- Al ° 



B\e~ A \ 
0, 



, g(y) 



/\6 ^>np 

The wave function appearing in (3.14) has the form 

V>i(A, x) ipi(X,x) 



e~ A ™B x 

r 







pn 

t„ 



pm 



-e- A \yc) 



*(A,x) = 



^ 2 (A, x) ip 2 (\x) 



where 



i>i (A, x) 
^2(A,x) 



as a; — > +oo and 

as x -> +oo 



-iAx r 



is the (m + n) x n Jost solution to (5.8) with the asymptotics 



V>i(A, x) 
^2 (A, x) 



is the (m + n) x to Jost solution to (5.8) with the asymptotics 



0, 
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Example 5.3 For the Zakharov-Shabat system in (5.1), the (right) Marchenko integral 
equations are given by 

'Q T (x + y)' 






M(x,y) + 



+ / dyM(x,z)Q T (z + y) 



y < x, 



M(x,y)- 



n r (x + y y 



dyM(x,z) Q r (z + y) 



t 



(5.12) 



y < x. 



The scalar potential u is recovered from the solution to (5.12) as 

0' 



1 



(5.13) 



(5.14) 



u(x) = 2[l Q]M(x,x) = -21(1,1)* 

By letting 

Q r (x + y)' 

a(x,y) = [M(x,y) M(x,y)], u(x,y) = 

_-n r (x + y)i 

we can transform (5.12) into (1.14), which is a 2 x 2 system of integral equations. Note 
that ui(x, y) given in (5.14) satisfies (1.10) with J = diag[l, —1]. As seen from (5.13), u is 
recovered from the solution to (1.14) as 

^0" 



u(x)=2[l 0]a(x,x) 



1 



In this case, f(x) and g(y) appearing in (1.13) are given by 

(\< J ' ' 1 " |f 

bU a * x o 



, g(y) 



e A ' y C} o 



o 



where A r is a constant p x p matrix with all eigenvalues having positive real parts, B r is 
a constant p x 1 matrix, and C r is a constant 1 x p matrix. In this case the wave function 
appearing in (3.14) is given by 



*(A,x) = 



(j>i(X,x) 4>i(X,x) 

02 (A, X) 02 (A, X) 



h(X,x) -0 2 (A*,x)* 
i 2 (A,x) 0i(A*,x)* 



where 
and 



01 (A, a; 

02 (A, X 

0i (A, x) 
2 (A,x) 



is the Jost solution to (5.1) with the asymptotics 



-iAx 







as x 



-oo 



is the solution behaving as 







+ o(l) as a; — > — oo. 
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Example 5.4 Example 5.3 is also valid for the Zakharov-Shabat system given in (5.8), 
where M, M, fi r , A r , £? r , C r , /, g have now sizes (m + n) x m, (m + n) x n, m x n, p x p, 
p x n, m x p, (m + n) x 2p, 2p x (m + n), respectively, and the eigenvalues of the constant 
matrix A r have all positive real parts, and B r and C r are constant matrices. The matrix 
system in (5.12) can be written as in (1.14), which is now an (m + n) x (m + n) system of 
integral equations The m x n potential matrix u is recovered from the solution to (1.14) 
as 

u(x) = 2 [I m mn ] a(x,x) 



0, 



In this case the wave function appearing in (3.14) is given by 



4>i{\x) 0i (A, x) 
02 (A, x) 4>2 (A, x) 



where 



0i (A, x) 
2 (A,x) 



as x — > — oo and 
0, 



is the (m+n) x m Jost solution to (5.8) with the asymptotics 

is the (m + n) x n Jost solution to (5.8) with the asymptotics 



-i\x 
Or 



6 Im 



ipi(X,x) 

V> 2 (A,£) 



J mn 



as x — > — oo. 



6 -in 

In the next three examples we discuss the Schrodinger equation on the full and half 
lines. In each example we have n bound states at k = ikj with the corresponding norming 
constant Cj. Here distinct positive numbers and ci , . . . , c n are positive. We 

begin with the familiar example of the Schrodinger equation on the full line [15] and then 
discuss the Schrodinger equation on the half line with various boundary conditions at x = 
[5]. 

Example 5.5 Consider the Schrodinger equation on the full line 

d 2 <p 



dx 2 



+ u(x) cp = k cp, x E (— oo, +oo), 



where u is a real- valued integrable potential with a finite first moment. The corresponding 
(left) Marchenko equation is a scalar integral equation and is given by (1.7), from whose 
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solution the potential u is recovered as 



u(x) = -2 



da(x, x) 
dx 



In this case, f(x) and g(y) appearing in (1.13) are given [5] by 



f(x) = [ae K1X ... c n e K ™ x ], g(y) = 



-Kiy 



-K n y 



where Cj is the norming constant for the bound state at k = iKj. The wave function 
appearing in (3.14) corresponds to the Jost solution from the left satisfying ^(k, x) = 
e lkx [l + o(l)] asi^ +oo. In this case (3.14) holds with the involution J being the scalar 
quantity equal to —1. 

Example 5.6 Consider the Schrodinger equation on the half line with the Dirichlet bound- 
ary condition at the origin, i.e. 



dx 2 



+ u(x) (p = k 2 tp, (p(k, 0) = 0. 



(5.15) 



The GePfand-Levitan integral equation arising in the related inverse scattering theory is 
given by 

a(x, y) + u(x, y) + / dya(x, z) u(z, y) = 0, < y < x, (5.16) 
Jo 

where the kernel u> is real valued and symmetric, i.e. lo(x, y) = u>(y, x), and thus it satisfies 
(1.10). Note that (5.16) is already in the form of (1.15). The potential is recovered as 



u(x) = 2 



da(x, x) 
dx 



In this case, f(x) and g(y) appearing in (1.13) are given by given [5] by 



f(x)= — sinh(Kix) ... — smh(K n x) , g(y) = 



— sinh(/tiy) 



sinh(K n y) 
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where Cj is the norming constant corresponding the bound state at k = ikj. In this case 
the wave function ty(k,x) is the regular solution to (5.15) satisfying the initial conditions 
^(k,0) = and \I/'(/c,0) = 1, where the prime denotes the cc-derivative. The relationship 
between the wave function ^(k, x) and a(x, y), instead of (3.14), is given by [5] 



x) 



sin(/cx) 
k 



+ / dya(x,y) 



sin(/cy) 
k 



< y < x, 



with the inverse transform given by 

1 f°° 

a{x,y) = — I dkk^(k,x) sin(/c,y), < y < x. 

K J-oo 

Example 5.7 Consider the Schrodinger equation on the half line with a selfadjoint bound- 
ary condition at the origin, i.e. 



dx'- 1 



+ u(x) if = k 2 (p, (p'(k, 0) + cot 9 ip(k, 0) = 0, 



(5.17) 



where 9 is a constant in the interval (0, n). The Gel'fand-Levitan integral equation arising 
in the related inverse scattering theory is given by 

r 

a(x,y) + u(x,y) + / dz a(x, z) u(z, y) = 0, < y < x, (5.18) 
Jo 

where the kernel u is real valued and symmetric, i.e. u>(x, y) = u{y, x), and thus it satisfies 
(1.10). Note that (5.18) is already in the form (1.15). The potential is recovered as 

da(x, x) 



u(x) 



dx 



In this case, f(x) and g(y) appearing in (1.13) are given [5] by 



ci cosh(fti2/) 



f(x) = [c 1 cosh(K 1 x) ... c n cosh(K n x) ] , g(y) 

cosh(n n y) 

where Cj is the norming constant for the bound state at k = iKj as in Example 5.6. In this 
case the wave function ^(k,x) is the regular solution to (5.17) satisfying the initial con- 
ditions \l/(/c,0) = 1 and ^'(k,0) = — cot 9. The relationship between ^(k : x) and a(x,y), 
instead of (3.14), is given by 



^(fc, x) = cos(kx) + J dy a(x, y) cos(/cy), < y < x, 
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with the inverse transform given by 

1 f°° 

ct(x, y) = — I dk ^(k, x) cos(ky), 0<y<x. 

7T ./-ex, 



6. DARBOUX TRANSFORM FOR THE ZAKHAROV-SHABAT SYSTEM 

In order to illustrate the significance of the results presented in this paper, we will 
derive the Darboux transformation for the Zakharov-Shabat system given in (5.1) when 
one bound state is added to the spectrum. In particular, to the spectrum of (5.1) we 
will add one bound state at A = Ai G C + with the norming constant ci, where c\ is a 
complex constant and we use C + to denote the upper half complex plane. The potential 
appearing in (5.1) will then change from u(x) to u(x) and the wave function appearing in 
(5.7) will change from #(A, x) to *(A, x). With the help of (3.20), (3.21), and (5.5), we will 
explicitly evaluate u{x) — u(x) and \I/(A, x) — \I/(A, x) in terms oiu(x), Ai, c±, and \I/(A, x). 
Note that u(x) and \I/(A, x) each consist of a one-parameter family of potentials and wave 
functions, respectively, with c\ being the parameter. From (5.7) we see that, in evaluating 

is 

as x — > +oo. At the end of this 



*(A, x) — \I/(A,;r), it is sufficient to obtain 
the Jost solution to (5.1) with the asymptotics 



ipi{\x) 




i>i{\x) 


, where 


ipi{\,x) 


_ip 2 {\x) _ 




_1p 2 {\,x) _ 


_1p 2 {\,x) _ 



iXx 



section we will compare our result with those in literature. 

Before we derive our Darboux transformation, we need some identities related to the 
solutions to (5.1). Let us use an overdot to indicate the derivative with respect to A. The 
following result is already known and its proof is omitted. Its proof can easily be obtained 
by using the integral representation of the Jost solution [30] . 





Proposition 6.1 Let 



ipi(X,x) 



be the Jost solution to (5.1) with the asymptotics 
as x — > +oo, where u is an integrable potential. If Ai G C + , then V'i(Ai,x), ^(Ai, x), 
■01 (Ai , x), and ip 2 (Ai, x) all vanish as x — > +oo. 

Even though the identities in the following proposition are known [1-3,30,35], we 
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provide a brief proof for the convenience of the reader. 
Proposition 6.2 Let 



be the Jost solution to (5.1) with the asymptotics 
as x — > +00, where u is an integrable potential. We then have the following identities 
ipi(X,x) -ip 2 (X,x) 







d_ 

dx 



d_ 

dx 



ip 2 (X,x) ipi(\,x)* 

ipi{X,x) ipi{X,x) 

1p 2 (X,x) tp2(\x) 

ipi(\,x) ipi(\i,x) 

1p2{\x) V2(Ai,x) 

ipi{X,x) -ip2{Xi,x)' 
ip 2 (X,x) i)i{\ux)* 
where Im[X] is used to denote the imaginary part of A 



d 

dx 



d 

dx 



2Im[A] [|^i(A,^)| 2 - |^ 2 (A,x)| 2 ] , (6.1) 

= 2iVi(A,x)^2(A,x), (6.2) 

-i(X - Ai)[V»i(A,x) ip2(Xi,x) + V'2(A,x)V'i(Ai,x)], (6.3) 

-i(A - A*) [^(A, x) Vi(Ai, x)* - V 2 (A, x)^ 2 (Ai, x)*], (6.4) 



PROOF: From (5.1) we obtain 

V4(A, x)* = zA*V>i(A, x)* + u(x)* x)*, 

ip' 2 (\, x)* = -zA*V> 2 (A, x)* - u(x) ipi(X, x)*, 

ip[(X, x) = -z^i(A, x) - iXipi(X, x) + u{x) V> 2 (A, x), 

ip' 2 (X, x) = iip 2 (X, x) + iXip 2 (X 7 x) - u{x)* ipi(X, x). 
The identities given in (6.1)-(6.4) are derived directly from (5.1), (6.5), and (6.6). | 

Using Propositions 6.1 and 6.2, we have the following result. 



(6.5) 



(6.6) 



Corollary 6.3 Let 



V>i(A,x) 
V> 2 (A,x) 



be the Jost solution to (5.1) with the asymptotics 



as x — > +00, where u is an integrable potential. We then have the following identities: 



dsip 1 (X 1 , s) V> 2 (Ai, s) = - ipi(\i,x) ip 2 (X!,x) - t/> 2 (Ai,x) ?/>i(Ai,x) 



|Vi(A 1 ,x)| 2 + |V^(Ai,x)| s 



poo 

£ d S (|^(A 1 , S )| 2 -|^ 2 (A 1 , S )| 2 )= 2Lii[Ai 

poo 

/ (is [V>2(Ai, s) V»i(A, s)+V»i(Ai, s) V> 2 (A, s)] 



(6.7) 
(6.8) 

(6.9) 



A- Ai 



[V>i(Ai,x) ij>2(\,x) - ^ 2 (Ai,x) V»i(A,x)] , 
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J X 



ds [V>i(Ai, s)* s) - ^ 2 (Ai, s)* ?/> 2 (A, s)] 



[V>i(Ai, x)* V»i(A, x) + ^(Ai, x)* ^2 (A, x) 



(6.10) 



A - A* 

Let us now use our systematic approach to derive the one-parameter family of the 
Darboux transformations for (5.1) when one bound state at A = Ai G C + with the norming 
constant c\ is added to the spectrum. Let us choose in (5.6) 

A\ = —iXi, B\ = 1, C\ = c\. 

Then, using (5.6) and (5.7) in (4.5), (4.7), and (4.10) we obtain 

n(x) = 



ciV>i(Ai,x) ^ 2 (Ai,x)* 
ci^ 2 (Ai,x) -ipi(\i,x)* 



q(x) 



^ 2 (Ai,x) Vi(Ai,x) 

c^Vi(Ai,x)* -c^2(Ai,x) s 

Ti(x) r 2 (x)' 

r 3 (x) r 4 (x) 



r(x) = 



where we have defined 



/•OO 

ri(x) :=T 4 (x)* := l + 2ci / ds V>i(Ai, s) V> 2 (Ai, s), 

J X 

r 3 (x) 



r 2 (x) := 



l c ll 



ds (|^ 1 (A 1 , S )| 2 -|V 2 (A 1 , S )| 2 ) 



Using (6.7) in (6.14) and using (6.8) in (6.15) we get 



r 3 (aO _ r , , _ |^i(Ai,x)| 2 + |^ 2 (A 1 ,x)| 2 



(6.11) 

(6.12) 
(6.13) 



(6.14) 
(6.15) 



r 4 (x)* = Ti(x) = 1 + %c\ ^i(Ai,x)^ 2 (Ai,x) -^ 2 (Ai,x)^i(Ai,x)j , (6.16) 



(6.17) 



Next, we present the main result in this section, namely, the Darboux transformation 
for the Zakharov-Shabat system given in (5.1). 
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Theorem 6.4 When one bound state at A = Ai G C + with the norming constant c\ is 
added to the spectrum of (5.1), the potential u(x) changes to u(x) and the Jost solution 



^2(A,x) 



changes to 



ijji (A, x) 
t/> 2 (A,x) 



u(x) — u(x) 



according to the following Darboux transformation 
P (x) 



ir^z)? + \ Cl \*r 2 (xy 



ipi(X,x) 




ipi(X,x) 


1 


Pi P2 




"P 5 " 


_^ 2 (A,x)_ 




_ip 2 (X,x)_ 


|ri(x)| 2 + | Cl | 2 r 2 (x) 2 


_Ps P 4 . 




. P 6. 



where we have Ti(x) andT 2 (x) are as in (6.16) and (6.17), respectively, and 
P :=2c 1 Vi(A 1 ,x) 2 r 1 (x) Ht -2 C n^ 2 (Ai,x)*] 2 r 1 (x) 
+ 4|ci|Vi(Ai,x)V2(Ai,x)*r 2 (x), 

Pi := -|c 1 | 2 V 2 (Ai,a;)*r 2 (a;)-c 1 ^ 1 (A 1 ,a;)r 1 (x)*, 
P 2 := |ci| 2 ip 1 (X 1 ,x) T 2 (x) -c* ip 2 (X 1 ,x)*T 1 (x), 
P3 := |ci| 2 ipi(\i,x)* r 2 (x) - Ci^ 2 (Ai, x) 
P4 := |ci| 2 ^ 2 (Ai,x) T 2 (x) + cj ^i(Ai,x)*ri(x), 



P 5 
P 6 := 



A- Ai 
— i 



[V>i(Ai,x) ij>2(\x) - ^ 2 (Ai, x) ^i(A,x)] , 



A - Aj 

PROOF: From (3.21) and (5.5) we have 



[V>i(Ai, x)* i) X (A, x) + ^ 2 (Ai, x)* -02 (A, x) ] 



w(x) — -u(x) = 2 [ 1 0]n(x)r(x) 1 q(x) 



(6.18) 



(6.19) 



(6.20) 



Thus, using (6.11)-(6.13), (6.16), and (6.17) in (6.20) we obtain (6.18). From (3.18), (3.20), 
(5.7), and (6.12) it follows that 



</>i(A,x) 




V>i(A,x) 


_^ 2 (A,x)_ 




_ip 2 {X,x)_ 



= —n(x) r(x) _1 / 

J X 



ds q(s) 



(6.21) 
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Using (6.11)-(6.13), (6.16), and (6.17) on the right hand side of (6.21), with the help of 
(6.9) and (6.10) we obtain (6.19). | 

Let us now compare our systematic approach to the Darboux transformation for the 
Zakharov-Shabat system given in Theorem 6.4 with some of the results in the literature. 
Some special cases of the Darboux transformations were obtained explicitly [10,23,24] for 
the Zakharov-Shabat system when a bound state is added at A = Ai, whereas our own 
method provides all the Darboux transformations in this case. When a bound state is 
added at A = Ai, one special Darboux transformation, named an elementary Darboux 
transform [10,23,24], corresponds to (cf. (1.9) of [10]) 



u(x) = = «'(*) + XXMx) - ^#4 «(*) 2 - 



(6.22) 



V»i(Ai,x)* " ' "" ^i(Ai,x) 

We will see that this Darboux transformation corresponds to the choice of 1/ci = in our 
one-parameter family given in (6.18). A second elementary Darboux transform [10,23,24] 
is given by (cf. (1.11) of [10]) 



u(x) = 4M = "«'(*) - 2i\{ <x) - ^jM u{xf 



(6.23) 



V»2(Ai,x) ' ' ^(Aijx)" 

which corresponds to an added bound state at A = A^. A combination of these two yields 
the following Darboux transformation (cf. (1.14) of [10]) when a bound state at A = Ai 
and another one at A = A^ are added to the spectrum of (5.1): 

4Im[Ai] ipi(Xi,x) ip2(Xi,x) 



u(x) - u(x) - , . 

|Vi(Ai,x)| 2 + |V 2 (Ai,x)| 2 

The Jost solution is then transformed according to (cf. (2.41) of [10]) 



(6.24) 



$i(X,x) 
$2(A,x) 



= | _t( \ — \ \ -l. Im[Ai]>g(Ai,x) 



V»i(A,x) 
ip 2 {\x) 



(6.25) 



where the matrix «S(Ai, x) is defined as 

IV^(Ai,x)| 2 



S(X u x) := 



-ipi(Xi,x) ^(Ai, x)' 
-ipi(Xi,x)*ip 2 (Xi,x) |V»i(Ai,x)| 2 
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Note that (6.25) implies that 



V>i(Ai,x) = ^ 2 (Ai,x) = 0. 

The norming constant c\ appearing in (6.18) is given by 

-1 



ci = 



i 

J — c 



(is V>i(Ai, s) t/> 2 (Ai,s) 



and hence the Darboux transformation given in (6.22) corresponds to a very particular 
choice of the norming constant ci, namely 1/ci = 0. Similarly, for the Darboux transfor- 
mation given in (6.23), the reciprocal of the norming constant is zero. Consequently, the 
Darboux transformation given in (6.24) corresponds to one particular choice among the 
two-parameter family u(x). 

We can also see that each of the three Darboux transforms given in (6.22)-(6.24) 
are very special by considering the easiest case where u(x) = 0. In that case (6.22) is 
understood as u{x) = and (6.23) and (6.24) each impose u(x) = 0. When u(x) = 0, the 



Jost solution 



^ 2 (A,x) 



to (5.1) with the asymptotics 

^i(A,x)=0, M\x)=e tXx 





iXx 



as x — > +oo is given by 

(6.26) 



Then, we see that (6.24) also imposes u(x) = 0. On the other hand, using (6.26) in (6.11), 
(6.12), (6.16), (6.17), and Theorem 6.4, our own procedure yields the one-parameter family 
of potentials and wave functions 

-8cJ (Im[Ai]) 2 e- 2iA i x 



u(x) = 
$i(A,x) = 



4(Im[Ai]) 2 + |ci| 2 e- 4Im [Ai]*' 
Aic\ (Im[Ai]) 2 e- 2lA i x+ax 



(6.27) 



(A- At) 4(Im[Ai]) 2 + |ci| 2 e- 4Im [Ai]x 



2i 


Cl 


2 (Im[Ai])e- 4Im [ Al ] x+iAl 


(A -At) 


4(Im[A!]) 2 + 


Cl 


2 e -4Im[Ai] x 
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Using [1-3,29] 

1 



+ o(l) 



T(A) 

where T(A) is the transmission coefficient corresponding to u(x), we get T(A) 



x — > — oo, 

A- 



A- Ai 

We note that the quantity given in (6.27) is the one-soliton potential and the corresponding 
reflection coefficients are zero. 

In [29] a Darboux transformation was studied for the Zakharov-Shabat system in 
two spatial variables, and by eliminating one of the variables, the corresponding Darboux 
transformation was provided in the form of an ordinary differential equation (cf. (4.9a) 
of [29]) involving u'(x), u'(x), u(x), and u(x). The same differential equation was derived 
earlier (cf. (4.7) of [26]) by Levi, Ragnisco, and Sym, who studied the equivalence of 
the dressing method [35] and the Darboux transformation for the Schrodinger equation 
and indicated that the result also holds for the Zakharov-Shabat system. These authors 
obtained a formula, which is given as (4.6) in their paper [26], connecting u(x) and u(x) 
through an intermediate function £ expressed in terms of the four entries of a matrix- 
valued solution to (5.1). A first-order differential equation similar to (4.9a) of [29] was 
earlier derived by Gerdzhikov and Kulish (cf. (14) of [19]). In [28], when N bound states 
are added to the spectrum, the change in the potential is expressed (cf (20) of [28]) in terms 
of ratios of the determinants of two 2N x 2N matrices differing only in their last columns; 
such matrices are constructed by determining the zeros of certain polynomial equations in 
A and by solving certain linear algebraic equations. Various other authors (see, e.g. [22]) 
presented similar formulas for the change in the potential when bound states are added to 
the spectrum. In most of these papers a "Darboux matrix" is constructed connecting the 
wave functions of the original and perturbed problems. 

One criticism of the result of [26] is that a matrix solution to (5.1) was evaluated (cf. 
(3.6) of [26]) at a A- value on the upper- half complex plane and also evaluated (cf. (3.8) of 
[26]) at a A- value on the lower- half complex plane. The same concern also applies to other 
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works (see, e.g. (19) of [28] and (3.8) of [22]). In general, we cannot expect the entries of a 
matrix solution to (5.1) to have extensions in A to both upper and lower complex planes, 
unless the class of potentials u(x) is very restrictive. 

A. APPENDIX: SOME NORM ESTIMATES 

In this appendix we derive some boundedness properties of the integral operators O 
defined in Section 1. 

Suppose u{x,y) is an N x N matrix function satisfying (1.2). Let 

poo poo 

Ni(x):=sup / dz\\uj(y,z)\\, N^x) := sup / dz \\oj(z, y)\\, 
y>x Jx y>x Jx 

so that (1.2) amounts to N 1 (x) + N oa (x) < +oo. 

Proposition A.l If the assumption in (1.2) is satisfied, then the integral operator O 
appearing in (1.3) is bounded on H^ xN with norm bound Ni(x) and on H^ xN with 
norm bound N OQ (x). 

PROOF: We directly verify that 

poo poo poo poo 

/ dy\\((Xl)(x,y)\\< dy dz \\P(x,z)\\ \\<j(z, y)\\ < N^x) / dz \\/3(x, z)\\, 

J X J x J x Jx 

and 

\\(pn)(x,y)\\ < JVoo^supH^y)!!, 

y>x 

which proves the proposition. | 

It is now clear that the (scalar) integral operator L u defined by 

poo 

(L u h)(y):= / dzh(z)\\oj(z,y)l 

J X 

is bounded on L 1 (x, +oo) with norm bound Ni(x) and on L°°(x, +oo) with norm bound 
N OQ (x). By the Riesz-Thorin interpolation theorem (cf. [36], Vol. II, Sec. XII. 1), L u is 
bounded on L v {x, +oo) for p e (1, +oo) with norm bounded above by 

N^xf^N^xf-^l 
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Since O maps a dense linear subspace of H^f xN (namely, its intersection with H^ xN H 
H™ xN ) into xN n H^ xN , with h(z) = \\/3(x, z)\\ we obtain the estimate 



l/p 

dz\\(pn)(x,z)\\*' 



< \\L u h\\ p < N^xf'PNcoixf-WriWhWp 

i/p 



^(xy/PN^x) 1 -^ 



dz\\(3(x,zW 



where || • || p is the L p -norm. Hence, is bounded on H^ IxN for p G (1, +00) as well. 
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